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We theoretically study effects of surface roughness on the magnetic response of small unconven¬ 
tional superconductors by solving the Eilenberger equation for the quassiclassical Green function 
and the Maxwell equation for the vector potential simultaneously and self-consistently. The param¬ 
agnetic phase of spin-singlet d-wave superconducting disks is drastically suppressed by the surface 
roughness, whereas that of spin-triplet p-wave disks is robust even in the presence of the roughness. 

Such difference derives from the orbital symmetry of paramagnetic odd-frequency Cooper pairs ap¬ 
pearing at the surface of disks. The orbital part of the paramagnetic pairing correlation is p-wave 
symmetry in the d-wave disks, whereas it is s-wave symmetry in the p-wave ones. Calculating the 
free-energy, we also confirm that the paramagnetic state is more stable than the normal state, which 
indicates a possibility of detecting the paramagnetic effect in experiments. Indeed our results are 
consistent with an experimental finding on high-Tc thin films. 

PACS numbers: 73.20.At, 73.20.Hb 


I. INTRODUCTION 

Diamagnetic response to an external magnetic field is 
a fundamental property of all superconductorsi. The 
Meissner current (coherent motion of the Cooper pairs) 
screens a weak magnetic field at a surface of a supercon¬ 
ductor. As a result, the phase coherence of superconduct¬ 
ing condensate is well preserved far away from the sur¬ 
face. A number of experiments, however, have reported 
the paramagnetic response of small superconductors and 
mesoscopic proximity structures 

Recent theoretical studies have suggested the existence 
of paramagnetic Cooper pairs in inhomogeneous super- 
conductors^^— . A spatial gradient of the superconduct¬ 
ing order parameter induces subdominant pairing corre¬ 
lations. The pairing symmetry of such induced Cooper 
pairs is different from that of principal Cooper pairs in 
bulk superconducting state — . For example, the prin¬ 
cipal Cooper pairs in high-U superconductors belong to 
the spin-singlet d-wave (even parity) class. In (110) di¬ 
rection of high-Tc cuprate, a surface acts as a pair breaker 
and suppresses the pair potential drastically. Simulta¬ 
neously, the spin-singlet odd-parity pairs are locally in¬ 
duced at the surface as a subdominant correlation. A sur¬ 
face generates odd-parity pairing correlations from the d - 
wave even-parity correlation because the surface breaks 
inversion symmetry locally. Since the pairing correla¬ 
tion function must be antisymmetric under the permu¬ 
tation of two electrons, the induced pairs have the odd- 
frequency symmetry—. To our knowledge, such induced 
odd-frequency pairs indicate the paramagnetic response 
to an external magnetic field. Odd-frequency Cooper 
pairs can be generated also from conventional supercon¬ 
ductors in the presence of spin-dependent potentialsii. 

In a previous paper—, we have shown that magnetic 
susceptibility of small enough unconventional supercon¬ 


ducting disks can be paramagnetic at a sufficiently low 
temperature. Odd-frequency Cooper pairs induced by 
a surface are responsible for the unusual paramagnetic 
Meissner effect. The magnetic response of Cooper pair is 
well characterized by so called “pair density” which is de¬ 
fined by diagonal elements of the response function to a 
magnetic field. Even-frequency Cooper pairs have a pos¬ 
itive pair density, whereas induced odd-frequency pairs 
have a negative pair density. So far an experiment has 
reported the decrease of pair density at low temperature 
in high-Tc superconducting films on which internal sur¬ 
faces are introduced by the heavy-ion irradiatioiii^. Thus 
our theoretical results are consistent with the experiment 
at least qualitatively. However, signs of the paramagnetic 
effect in the experiment — is much weaker than our theo¬ 
retical prediction. The discrepancy may come from sam¬ 
ple quality at surfaces. Artificially introduced internal 
surfaces can be very rough in the experiment, whereas 
surfaces are specular in the theory. Actually, several the¬ 
ories have pointed out that the surface roughness affects 
properties of the surface Andreev bound states of a high- 
Tc superconductor 

The purpose of this paper is to clarify effects of sur¬ 
face roughness on the paramagnetic Meissner response 
of small unconventional superconductors. We consider a 
two-dimensional superconducting disk with spin-singlet 
d-wave symmetry or spin-triplet p-wave one. In numeri¬ 
cal simulation, we solve the Eilenberger equation and the 
Maxwell equation simultaneously and self-consistently. 
Surface roughness is considered through an impurity self¬ 
energy within the self-consistent Born approximation. 
We find that the surface roughness suppresses drastically 
the paramagnetic response of a spin-singlet d-wave su¬ 
perconducting disk. On the other hand in a spin-triplet 
p-wave disk, the paramagnetic property is robust even 
in the presence of surface roughness. The induced odd- 
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FIG. 1. (a) Schematic figure of a superconducting disk with 

rough surface, with R and w are the radius of the disk and 
the width of the disordered region, respectively. An external 
magnetic field is applied in the z direction. The origin of a 
spatial coordinate is located at the center of the disk. The 
pair potential in momentum space for the d-wave supercon¬ 
ductor and that for p-wave one are illustrated in (b) and (c), 
respectively. 


frequency pairing correlation has p-wave symmetry in 
the former, whereas it has s - wave symmetry in the lat¬ 
ter. In addition, we also confirm that the paramagnetic 
superconducting states are more stable than the normal 
state by calculating free-energies. 

This paper is organized as follows. In Sec. II, we ex¬ 
plain the theoretical method to analyze the magnetic re¬ 
sponse of small superconducting disks. In Sec. Ill, we 
discuss the magnetic response of small superconducting 
disks with a rough surface. In Sec. IV, we discuss the 
stability of a paramagnetic state by calculating the free- 
energies of a superconducting state. We summarize this 
paper in Sec. V. 


II. FORMULATION 


Let us consider a superconducting disk in two- 
dimension as shown in Fig.[TJa), where R is the radius of 
the disk. To describe rough surfaces, we introduce ran¬ 
dom impurity potentials near the surface. The width of 
the disordered region w is measured from the surface as 
shown in Fig. [ija). An external magnetic field is 

applied in the z direction. Throughout this paper, we 
use a unit of h = c = ks = 1 with ks and c being the 
Boltzmann constant and the speed of light, respectively. 

Superconducting states in equilibrium are described by 
solutions of the quasiclassical Eilenberger equation^^, 

ivpk ■ Vr g-k -h t:, g] = 0, (1) 


where uj’ is the Fermi velocity, k is the unit vector on 
the Fermi surface, g and H are defined as follows, 
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where = (2n -|- l)7rT is the fermionic Matsubara fre¬ 
quency, n is an integer number, and T is a temperature. 
In this paper, the symbol • ”• represents a 4 x 4 matrix 
structure, •" • represents a 2 x 2 matrix structure in spin 
space and tro is the identity matrix in spin space. A 
vector potential is denoted by A{r). We introduced a 
definition K{r, k, iujn) = K*{r, —k, iujn) for all functions 
iL(r, k, iuJn)- Effects of rough surfaces are taken into ac¬ 
count through an impurity self-energy of a quasiparticle 
defined by, 

i f dk 

t{r,iuJn) = Q{\r\-R-\-w) — / — g{r,k,iujn), 

( 5 ) 

where tq is the life time of a quasiparticle and 0(a;) is the 
Heviside step function. The mean free path of a quasi¬ 
particle is defined by £ = vpTo in the disordered region. 
The anomalous Green function f{r,k,iujn) is originally 
defined by an average of two annihilation operators of an 
electron. The relation 

f{r,k,iuJn) =-f'^{r,-k,-iuJn), ( 6 ) 

represents the antisymmetric property of the anomalous 
Green function under the permutation of two electrons, 
where T represents the transpose of a matrix. 

The direction of k in two-dimensional momentum 
space is represented by an angle 6 measured from the 
X axis, (i.e., kx = cos 9 and ky = sin0). In what fol¬ 
lows, we consider two unconventional superconductors 
with different pairing symmetries. One is spin-singlet d- 
wave symmetry A{r,9) = A(r) sin(20) Z(T 2 . The other is 
spin-triplet p-wave symmetry A{r,9) = A(r) cos(0) di, 
where dj for j = 1-3 are the Pauli matrices in spin space. 
A d-wave and a p-wave pair potentials in momentum 
space are shown schematically in Fig. [TJb) and (c), re¬ 
spectively. We do not consider any spin-dependent po¬ 
tentials in this paper. The matrix structure of Green 
function is represented by 


g{r,e,iujn) =g{r,9,iuJn)do, 


( 7 ) 


f{r,e,iujn) =f{r,e,iujn) x 


(T 2 for spin-singlet 
—idi for spin-triplet ’ 
( 8 ) 


with scalar Green functions g{r,9,iuin) and f{r,9,iujn)- 
Spatial dependence of A(r) is determined self- 
consistently from the gap equation 

/ dG 

— f(r,e,iu}n)Vx{9), (9) 


where Nq is the density of states per spin of normal metal 
at the Fermi level, go is the coupling constant, and 14 
represents attractive electron-electron interactions with 
X = p-wave or d-wave indicating pairing symmetries. 
The interaction kernel 14 depends upon pairing symme¬ 
tries as 


Vx{9) = 


^{r, fc, iuJn) = {ioJn + evpk ■ A{r)) do, 


( 4 ) 


2 cos 6 for X = p-wave 
2sin(20) for x = d-wave' 


( 10 ) 
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FIG. 2. Local susceptibilities of the small superconducting disks. The results of a d-wave and those of a p-wave snpercondnctor 
with a clean surface (i.e., w = 0, = 0) are presented in (a) and (b), respectively. The results of a d-wave and those of a 

p-wave superconductor with rough surface (w = 3^o, = 1-0) are demonstrated in (c) and (d), respectively. The parameters 

used in the simulation are R = lO^o, Al = 5^o, tuc = lOAo, Lf®** = O.OlIdci and T = O.lTc. 


The constant iVogo is determined by 


{Nog) ^ = In 


+ 


E 


0<n<u;c/27rT 


.+ 1/2 ’ 


( 11 ) 


with Tc and uic being the transition temperature and the 
cut-off energy, respectively. 

In a type-II superconductor, an electric current is rep¬ 
resented by 


3{r) 


neypNo 

2i 


^E/£Tr[T3fcg(r,0, (12) 


with T 3 = diag[(T 0 ) —do]. From Eq. (TT^ and the Maxwell 
equation VxFf(r) = 47 rj(r), we obtain spatial profiles of 
a vector potential A{r) and a local magnetic field H{r). 
A local magnetic susceptibility is defined by 


1 H{r) - 
Hr ’ 


(13) 


where is an amplitude of external magnetic field 

applied in the z direction. By integrating the local sus¬ 
ceptibility, we obtain a susceptibility of a whole disk X 
as 


A = 


1 


'kl<« 


drx{r). 


(14) 


To solve the Eilenberger equation Eq. o in a disk 
geometry, we use a Riccati parametrization^— and 
a numerical method discussed in Ref. [ 2 ^. Using the 
parametrization, the Eilenberger equation can be sepa¬ 
rated into two Riccati type differential equations. When 
we solve the Riccati type equation along a long enough 
quasiclassical trajectory, solutions of the equation do not 
depend on initial conditional^. In this paper, the length 
of classical trajectories is more than 30 times of the co¬ 
herence length. Solving the Eilenberger equation and 
the Maxwell equation, we determine the pair potential, 
the vector potential, and the self-energy self-consistently 
with one another. At surfaces, we consider specular clas¬ 
sical trajectories for calculating the Green functions^S. 
The vector potential out side of the superconducting disk 
is A{r) = /2){—yx + xy) which gives a uniform 

magnetic field in the z direction, where x and y are the 
unit vector in the x and y direction, respectively. 


III. RESULTS 

Throughout this paper, we fix several parameters as 
R = lO^o and oJc = IOAq, where Aq is the ampli¬ 
tude of the pair potential at the zero temperature, and 
^0 = hvF/2,TrTc is the coherence length. Strength of the 
disorder is tuned by changing a parameter ^q/£. Width 
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FIG. 3. Local susceptibilities on the x axis (i.e., y = 0). The 
results for a d -wave disk and for a p -wave disk are plotted in 
(a) and (b), respectively. The all parameters are set to be the 
same as those of Fig.O The solid and the broken lines indicate 
the result for the disk with clean surface and those with the 
rough surface, respectively. The shadowed areas indicate the 
disordered regions. 
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FIG. 4. Amplitudes of pair potentials in real space for several 
azimuthal angle a measured from the x axis. The results of 
a d-wave disk with £/^o = 0.0 and £/^o = 1.0 are shown in 
(a) and (c), respectively. Since the pair potential is four-fold 
symmetric, we plot its spatial profile for a = 0, vr/S, and 7r/4. 
The results of a p-wave disk with £/^o = 0.0 and £/^o — 1.0 
are shown in (b) and (d), respectively. The pair potential for 
a — 0, tt/S, 7r/4, Stt/S, 7r/2 (from bottom to top) are shown 
under the two-fold symmetry of results. The pair potentials 
are normalized to ABuik(r). The all parameters are set to be 
the same as those of Fig. [2l The shadowed areas in (c) and 
(d) indicate the disordered regions. 


of disordered region is w = 3^o because odd-frequency 
pairs induced by a surface localize within this range. All 
lengths are measured in units of ^o- The current density 
is normalized to Jq = ?ic^/ 47 r|e|^o. Here we express h 
and c explicitly to avoid misunderstandings. The charac¬ 
teristic length scale of the Maxwell equation is the pen¬ 
etration depth defined as Ap = and 

is a parameter in numerical simulations. In this paper, 
we choose Xl = 5^o to realize type-II superconductors 
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FIG. 5. Spatial dependencies of the decomposed pairing 
functions at the lowest Matsubara frequency on the a;-axis for 
(a) a clean d-wave, (b) a clean p-wave, (c) a disordered d- 
wave, and (d) a disordered p -wave disk. The parameters are 
set to be the same values as those in Fig. [2l The shadowed 
areas indicate the disordered regions. 


and fix = O.OIiJca- Here iJca = ^c/|e|^Q is the 

second critical magnetic field. The first critical mag¬ 
netic field at low temperature is estimated as Hci = 
HcA^o/XL?log{XL/^o) - 0.03He2 > at Ap = 5^0- 
Thus vortices are not expected at low temperature. We 
start numerical simulations with initial conditions of a 
spatially uniform pair potential A(r) = |ABuik(T)| and a 
homogeneous magnetic field A(r) = /2){—yx+xy), 

where ABuik(T) is the pair potential obtained in a ho¬ 
mogeneous superconductor at a temperature T. If we 
choose an alternative initial condition hosting a vortex 
in a superconductor, a vortex state might be realized in 
numerical simulations^ even for < Hd- Such vor¬ 
tex issue, however, goes beyond the scope of this paper. 

First we discuss calculated results for a small super¬ 
conducting disk with a specular surface (i.e., w = 0 , 
^o/£ = 0). Then we discuss effects of surface roughness 
by comparing numerical results in a disk with a rough 
surface with those with a specular one. 


A. Disks with a specular surface 

We briefly explain the paramagnetic Meissner effect 
in a superconducting disk with a specular surface. Fig¬ 
ure [ 2 Ka) shows the local susceptibility in a d-wave super¬ 
conductor. The result are four-fold symmetric reflecting 
the d -wave pair potential. The magnetic susceptibility is 
positive (paramagnetic) near four surfaces in the x and 
the y directions. In Fig. EJa), we show a spatial profile 
of the local susceptibility of Fig. [2ja) along the x axis 
at y = 0. We also show the amplitudes of pair poten¬ 
tial in real space |A(r)| in Fig. HDa). The pair poten¬ 
tial is calculated along a trajectory pa oriented by an 
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angle a measured from the x axis as shown iir Fig. |4l 
In a d-wave disk, the pair potential is four-fold sym¬ 
metric. The pair potential is strongly suppressed at the 
four surfaces in x aird y directioirs as a result of appear¬ 
ing surface Andreev bound states^, whereas it is totally 
constant in the directions in which the surface Andreev 
bound states are absent. In coirtrast to the pair poten¬ 
tial around a vortex core^, the results in Fig. SJa) are 
anisotropic in real space. As shown in Ref. [l^ the spatial 
gradient of pair potential generates paramagnetic pairing 
correlations with odd-frequency symmetry. Thus odd- 
frequency pairs are expected at four surfaces in the x and 
y directions, which explains inhomogeneous and angular 
anisotropic paramagnetic response in Fig. [2ja). In this 
paper, we analyze the frequency symmetries of Cooper 
pairs by decomposing pairing functions into a series of 
Fourier components. In a d-wave disk, the anomalous 
Green function are described by two components 

f{r, 9, iion) = [/ep(r', 0, iiOn) + /op(r', 9, *w„)] <72, (15) 

where /ep is an even-parity (d -wave) function represent¬ 
ing the principal pairing correlation and /op is an odd- 
parity function representing the induced pairing compo¬ 
nent at a surface. To satisfy Eq. (0), /op must be an odd 


function of ujri- We decompose pair functions f{r. 
as 

1 0 , iUJri) 


Vi{r,iLOn) = 2^Cf + Sl 

(16) 


Si= J ^ Re [f{r, 9, iujn)] sin (19), 

(17) 


f rIO 

Ci=J — Re [f{r, 9, iujn)] cos {19), 

(18) 


where I = 0,1,2 and 3 correspond to s-, p-, d- and /- 
wave orbital functions, respectively. In the presence of 
a magnetic field, the imaginary part of /(r, 9, icun) is in¬ 
duced by the vector potential as analytically shown in 
Appendix. We focus only on the real part of / to ana¬ 
lyze pairing symmetries. Figure [Sja) indicate the spatial 
profile of Vi{x, iuJo) at the lowest Matsubara frequency as 
a function of a; at y = 0. The pairing functions of induced 
Cooper pairs have p - and / -wave symmetries and their 
amplitudes localize near the surface. Such odd-frequency 
Cooper pairs shows the paramagnetic response to a mag¬ 
netic field. The surface also generates spin-singlet s-wave 
correlation. Its amplitude, however, is too small to con¬ 
firm at a scale of plot in Fig.[5l)a). Figure [6l)a) shows the 
spatial distribution of electric current on a d -wave disk. 
The diamagnetic current flows at the central region be¬ 
cause of the usual Meissner effect. Near the surfaces in 
the X and y directions, however, the current flows the 
opposite direction to the Meissner current. Therefore a 
small d -wave superconductor can be paramagnetic due to 
the induced odd-frequency Cooper pairs at its surface. 

A spin-triplet p-wave disk also indicates the similar 
paramagnetic effect as shown the results of magnetic sus¬ 
ceptibility in Figs. [2Kb), its spatial profile on the x axis in 


(a) clean d-wave 


(b) clean p -wave 
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(c) disordered d -wave (d) disordered p -wave 



x/(o x/io 


FIG. 6. Spatial distribution of current density on a super¬ 
conducting disk. The results of a d-wave and a p-wave su¬ 
perconductor with clean surface {w = Q, 1 = 0.0) are shown 
in (a) and (b), respectively. The results with rough surface 
are shown in (c) and (d). The outer circles indicate the edge 
of the disk and the inner circles in (c) and (d) indicate the 
effective interface between the clean region and the disordered 
region. The parameters are set to be the same values as those 
in Fig. [2] 


ISb), pair potential injJKb) and electric current in|5Kb). 
The results are two-fold symmetric reflecting the p -wave 
superconducting pair potential. The paramagnetic effect 
can be seen near the surface in the x direction because 
of induced odd-frequency Cooper pairs. The anomalous 
Green function is represented by Eq. (ITSl) with replacing 
d '2 by —idi- In the p-wave case, /ep represents induced 
pairing correlations and is an odd function of a;„. As 
shown in Fig.jSKb), /ep mainly consists of s- and d-wave 
pairing correlations. 

In the end of this subsection, we summarize an im¬ 
portant difference between the paramagnetic effect of a 
d-wave superconductor and that of a p-wave one. In 
a d-wave disk, surface odd-frequency Cooper pairs have 
a p- or /-wave symmetr y^^d^ ^ ^ p-wave disk, on the 
other hand, s - or d -wave odd-frequency Cooper pairs are 
responsible for the paramagnetic effec t^^d^ . In the next 
subsection, we will show that the paramagnetic response 
of a disk with rough surface depends sensitively on the 
orbital symmetry of induced odd-frequency pairs at the 
surface. 
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B. Disks with a rough surface 

Next, we discuss effects of the surface roughness on the 
magnetic response of a small superconductor. The calcu¬ 
lated results of the local susceptibility for a d -wave super¬ 
conducting disk with rough surface are shown in Fig. He), 
where we choose = l-O- Comparing Fig. [2])a) with 
(c), we can find that the surface roughness completely 
suppresses the paramagnetic response at the four sur¬ 
faces in a d-wave disk. The central region of the disk 
with the rough surface recovers the usual diamagnetic re¬ 
sponse. Such effect is demonstrated more clearly in a spa¬ 
tial profile of local susceptibilities at y = 0 in Fig. H®-)) 
where the shadowed area indicates the disordered region. 
The amplitude of pair potential in real space is shown in 
Fig. He). The pair potential in the disordered region is 
totally suppressed because the random impurity poten¬ 
tial acts as a pair breaker for d-wave Cooper pairs. Spa¬ 
tial profiles of decomposed pairing functions are shown 
in Fig. H^)- In the disordered region, a d-wave pair¬ 
ing function Vd is drastically suppressed due to impurity 
scatterings. The disordered region can be considered as 
a diffusive normal metal because the spatial profile of 
order parameter is proportional to Vd- Odd-frequency 
Cooper pairs are also fragile in the presence of surface 
roughness because they have p- or f -wave pairing sym¬ 
metry. Therefore, both the paramagnetic current and 
the diamagnetic one disappear in the disordered region 
as shown in Fig. He). The magnetic property of a disk 
is determined by that at the clean central region where 
even-frequency d -wave Cooper pairs stay and contribute 
to the diamagnetic response. We conclude that the para¬ 
magnetic effect in a d -wave disk is fragile in the presence 
of surface roughness because odd-frequency pairs have a 
p -wave or a / -wave orbital symmetry. 

A p-wave disk indicates qualitatively different mag¬ 
netic response from a d-wave one. The local suscepti¬ 
bility of a p-wave disk with rough surface is shown in 
Fig. Hd) and Fig. Hb) with a broken line. Although 
the surface is rough enough, a p-wave superconducting 
disk still show the strong paramagnetic response. The 
peak of the x in Fig- Hb) shifts to inside of the disk 
in the presence of surface roughness. This suggests that 
the Andreev bound states appear at a boundary between 
the clean central region and the disordered surface re¬ 
gion. Such Andreev bound states always accompany the 
paramagnetic odd-frequency Cooper pairs. In addition to 
this, the paramagnetic response in Fig.^b) suggests the 
penetration of odd-frequency Cooper pairs into the sur¬ 
face disordered region. The spatial profiles of the electric 
current in Fig.[6])d) shows that the paramagnetic current 
flows not only in the clean region but also in the disor¬ 
dered one. Odd-frequency pairs in a p-wave disk survive 
even in the presence of surface roughness because they 
have s -wave orbital symmetry. Therefore, the paramag¬ 
netic effect in a p -wave disk is robust against the surface 
roughness. 



T/T, T/T, 

FIG. 7. Temperature dependencies of the whole disk suscep¬ 
tibility X(T) for the d-wave disk (a) and the p-wave disk 
(b). 


C. Temperature dependence 

Here we discuss the magnetic susceptibility of a whole 
disk which is a measurable value in experiments. The 
disk susceptibility in a d-wave superconductor and that 
in a p -wave one are plotted as a function of temperature 
in Fig. Ha) and (b), respectively. We present the results 
for several choices of the disorder /^. In simulation, we 
first calculate the pair potential and the vector potential 
self-consistently at a temperature just below Tc under an 
external magnetic field = 0.01idc2. Then the tem¬ 
perature is decreased with keeping unchanged. In 
the clean limit = 0.0), both a d-wave disk and a p- 
wave one show the usual diamagnetic response just below 
Tc- With decreasing the temperature, however, the sign 
of susceptibility changes around T = Tp 0.3Tc for both 
cases. Here, we define Tp as a diamagnetic-paramagnetic 
crossover temperature. Blow Tp, superconducting disks 
show the anomalous paramagnetic response. The para¬ 
magnetic effect is stronger in lower temperature because 
odd-frequency Cooper pairs energetically localize at the 
Fermi level. 

In a d-wave disk, the reentrance is slightly suppressed 
in the presence of the moderate surface roughness with 
^o/£ = 0.1 as shown in Fig. H^')- When we increase 
the degree of the roughness further, the paramagnetic re¬ 
sponse gradually becomes weaker. At ^o/£ = 0.5, the re¬ 
sponse is diamagnetic and the susceptibility recovers the 
monotonic temperature dependence which is usually ob¬ 
served in large enough superconductors in experiments. 
In the d -wave superconductors, the specular Andreev re¬ 
flection is necessary for forming the surface bound states 
at the zero-energy^ and for appearing odd-frequency 
pairs. In other wards, odd-frequency pairs have p- or 
/-wave orbital symmetry. Therefore the rough surface 
brakes odd-frequency pairs and suppresses the paramag- 
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netic response. This conclusion is totally consistent with 
the experiments^, where temperature dependences of the 
pair density on a high-Tc superconducting film show a 
small reentrant behavior at low temperature. But the 
total pair density remains positive. Actually, the experi¬ 
mental data are very similar to the results for ^o/i = 0.2 
in Fig.[71[a). The experimental results can be interpreted 
as an appearance of a small amount of odd-frequency 
pairs. In the experiment, the surface roughness may par¬ 
tially breaks odd-frequency pairs because a number of 
the internal surfaces are introduced by the heavy ion ir¬ 
radiation. 

On the contrary to a d-wave disks, the susceptibil¬ 
ity of a p-wave disk X{T) shows the reentrance and 
the crossover to the paramagnetic phase at low temper¬ 
ature for all ^o/i as shown in Fig. [71[b). It has been 
pointed out that the surface Andreev bound states of a 
p -wave superconductor are robust under potential disor¬ 
der because of the pure chiral property of surface bound 
states^. In other wards, odd-frequency pairs accompa¬ 
nied by the Andreev bound states have s-wave orbital 
symmetryi^ii^. Since s-wave pairs is robust under the 
disordered potential, the paramagnetic effect of the p- 
wave superconductors persists even in the presence of 
surface roughness. We conclude that the robust param¬ 
agnetic response in a small size sample is a unique prop¬ 
erty of spin-triplet p-wave superconductors. Such prop¬ 
erty would enable us to identify the spin-triplet p-wave 
superconductivity in experiments. 


IV. STABILITY OF PARAMAGNETIC 
SUPERCONDUCTING STATES 


Generally speaking, a superconducting phase is more 
stable than a normal one as far as a superconductor is dia¬ 
magnetic and homogeneous^. Therefore a homogeneous 
paramagnetic superconducting phase is usually unstable. 
The calculated results in Sec. Ill, however, show that 
the paramagnetic phase on a small superconducting disk 
is spatially inhomogeneous. In such situation, it would 
be worthy to check if the paramagnetic phase is a sta¬ 
ble state at a free-energy minimum or a metastable state 
corresponding to a free-energy local minimum. In this 
section, we discuss the stability of paramagnetic phase 
in small unconventional superconductors by calculating 
the free-energy in clean superconducting disks. 

The free-energy is calculated from the quasiclassical 


Green functions^, 


Fs - 

Fn = J dr T{r), 

(19) 

Hr) 

= •Fa(t’) + Tnir), 

( 20 ) 

Tnir 

^ {H{r) - 

8tt 

( 21 ) 

THr 

)=Tfir)+Tg{r), 

( 22 ) 

Hir] 

\=ttNo f ^T^X*{r,e)f{r,e,iuJn), 

(23) 




X / duj Re{g{r,9,iuj) — 1}, (24) 

J LOn 

where is the condensation energy density of elec¬ 

tron system and Tnif') is the energy density of a mag¬ 
netic field. We introduce an additional energy cut-off 
ujc 2 to evaluate the integration in Eq. (IMl) . In this pa¬ 
per, we set ujc 2 = 400Ao so that f drT[r) reaches to a 
converged value. The free-energy densities normalized to 
To = Vo|Aop /2 which is the condensation energy den¬ 
sity in a homogeneous s-wave superconductor. A tem¬ 
perature is set to be sufficiently low at T = O.lTc so that 
a superconductor is in the paramagnetic phase. 

The calculated results of the free-energy density for a 
d-wave disk and those for a p-wave one are shown in 
Fig. UKa) and (b), respectively. The free-energy density 
is calculated along a trajectory pa oriented by an angle a 
which is measured from the x axis as shown in Fig. [ 8 ] In 
a d-wave disk, the results in Fig. [5Da) show that T{r) is 
negative around the disk center. However it becomes pos¬ 
itive near the surfaces at a = 0. On the other hand, T{r) 
at a = 7 r /4 is almost flat and is always negative along 
the trajectory because odd-frequency pairs are absent in 
this direction. The results for a = 7 r /2 are identical to 
those for a = 0 due to the four-fold symmetry. The 
free-energy density varies gradually from the line with 
a = 0 to that with a = 7 r /4 as increasing a from a = 0. 
The free-energy of whole disk J drT{r) can be negative 
even though the disk is in the paramagnetic phase. In 
a p-wave disk, the increasing of T{r) occur only in two 
surfaces due to its p-wave symmetry. Therefore, we can 
conclude that both d- and p-wave superconducting states 
are more stable than a normal state even in their para¬ 
magnetic phase. 

To analyse the details of the free-energy density fur¬ 
ther, we decompose the free-energy density at a = 0 
into Th and J'a as shown in Fig. [IKa). When a super¬ 
conductor shows perfect Meissner effect, the energy of a 
magnetic field becomes Th = (id®’'*)^/ 87 r. In a small su¬ 
perconductor (i.e., R ^ Al), an external magnetic field 
penetrates into whole the disk. This suppresses Th from 
(p^ext) 2 / 87 r at the center of the disk. Near the surface, 
on the other hand, paramagnetic odd-frequency Gooper 
pairs attract a magnetic field, which increases Th locally. 




Pa/Co 



FIG. 8. Free-energy densities on a trajectory pa of a d-wave 
disk (a) and a p-wave disk (b), where T{r) is normalized to 
To — A^oAo/2. The angle a is measured from the x axis 
as shown in schematics in the figures. The temperature is 
set to T — O.lTc so that the superconducting disks show the 
paramagnetic response. The second energy cut-off is set to 
ljJc2 = 400 Aq. The other parameters are fixed as R = lO^o, 
Ai = 5Co, = 10Ao, = O.OlHc,. 


The appearance of odd-frequency pairs also increases Ta 
as discussed in Appendix. As a result, T{r) becomes pos¬ 
itive at the disk surface in both Fig. |9ja) and (b). The 
condensation energy Ta is negative at the disk center. 
While, near the surface, Ta increases due to the sup¬ 
pression of the pair potential there in both Fig. EJa) and 
(b). In appendix, we analytically calculate the Green 
functions and the free-energy density near the surface 
of a semi-infinite px-'Vf&Ye superconductor. The results 
in Eq. (IA.21I) indicates that the free-energy density is 
positive due to the appearance of odd-frequency Cooper 
pairs. The numerical results in Fig. m a) and (b) show 
that the free-energy of whole disk f dr J^(r) remains neg¬ 
ative because odd-frequency pairs are confined only near 
the surface. Therefore the paramagnetic superconduct¬ 
ing state on d - and p -wave disks is more stable than the 
normal state. 

Next, we study energetic properties of odd-frequency 
pairs. In our simulation, it is possible to obtain two 
superconducting states: a superconducting state in the 
absence of a magnetic field = 0 and that in the 
presence of a magnetic field ^ 0. At = 0, 

there is no electric current everywhere. Superconduct¬ 
ing states at ^ 0, on the other hand, carry electric 



x /^0 





FIG. 9. Condensation energy density and energy density of 
a magnetic field for a = 0. The results of a d-wave disk and 
those of a p-wave disk are plotted in (a) and (b), respectively. 
T, Ta, and Th are indicated in solid lines, broken lines, and 
dotted lines, respectively. All of the energy densities are nor¬ 
malized to To = NoAo/2, which is the condensation energy 
density of a homogeneous s-wave superconductor. The all 
parameters are set to be the same as those of Fig. |8] 


currents as shown in Fig. 5(a)-(d). Here we compare con¬ 
densation energies of such two different superconducting 
states as shown in Fig. [TOl where we plot Ta on a p- 
wave disk for a = 0. The solid line and the broken line 
indicate results for ^ 0 and those for iJ®’'* = 0, 

respectively. As shown in Fig. 1101 Ta near surfaces for 
jjext ^ Q jg lower than that for = 0. Odd-frequency 
pairing state in the presence of electric currents is more 
stable than that in the absence of electric currents. This 
energetic property explains the paramagnetic property 
of odd-frequency Cooper pairs. The argument above is 
valid also for a d-wave disk. In Appendix, we present 
analytical expression of difference between free-energy at 
ppext _ Q ^ results show that 

a magnetic field decreases the free-energy at low tem¬ 
perature because odd-frequency Cooper pairs have the 
paramagnetic property. 

V. CONCLUSION 

We have theoretically studied effects of surface rough¬ 
ness on the anomalous paramagnetic response of small 
unconventional superconducting disks by using the qua- 
siclassical Green function method. We conclude that the 
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FIG. 10. Condensation energy density Ta obtained with a 
external field (solid line) and without a external field (dotted 
line). Calculations are carried out for a p-wave superconduct¬ 
ing disk. The all parameters are set to be the same as those 
of Fig. H 


paramagnetic property of p -wave superconductors is ro¬ 
bust under the surface roughness because the p-wave 
superconductors host the s -wave odd-frequency Cooper 
pairs at their surface. On the other hand, the param¬ 
agnetic property in d-wave superconductor is fragile in 
the presence of the surface roughness. In this case, the 
odd-frequency pairs at the surface have p-wave orbital 
symmetry. We have also confirmed that the paramag¬ 
netic superconducting phase are more stable than the 
normal state by calculating the free-energy. 
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Appendix: Analysis in a semi-infinite pa;-wave 
superconductor 

In a semi-infinite superconductor in two-dimension, it 
is possible to obtain analytic expression of Green func¬ 
tions in the clean limit. The evaluation of a electric 
current and free-energy by using analytical expressions 
would be helpful to understand numerical results in the 
text. 

We assume that a superconductor occupies x > 0 and 
uniform in the y direction. An magnetic field applied 


in the z direction and its vector potential is given by 
A = A{x)y. The Eilenberger equation in 2 x 2 Nambu 
space reads. 


ivpk ■ Vg -I- [H,g] = 0, 


g{x,k,iujn) 


9 f 
spf -g 




H = 


iujn + evpk-A fA(x, fe) 

iSpA(x,k) —iujn — evpk 



( 1 even-parity 
1—1 odd-parity 


(A.l) 

(A.2) 

(A.3) 

(A.4) 


where a factor Sp depends on the parity of order param¬ 
eter and the Green functions satisfy g^ + Spff = 1. The 
Green functions can be expanded with respect to the vec¬ 
tor potential as 

9 =9^°^ + {-ievpk ■ A)du,„g^°'> 

+ ^{-ievpk ■ Afdl^g^°'> H-, (A.5) 

f + i-ievpk ■ 

+ ^{-ievpk ■ H-, (A.6) 

because a vector potential shifts the Matsubara fre¬ 
quency, where g^^^ and are the Green function in 
the absence of a vector potential. In what follows, we 
omit “(0)” from the Green function for simplicity. We 
note in Eq. (jA.6() that parity and frequency symmetry of 
the second term on the right-hand side are opposite to 
those of the first term because fc is a odd-parity function 
and dai^f changes a frequency symmetry!^. The imagi¬ 
nary part of an anomalous Green function represents a 
pairing correlation deformed by a vector potential. 

In the case of a p^-wave superconductor, it is possible 
to obtain a reasonable solution of the Eilenberger equa¬ 
tion at A = 0. When we assume the spatial dependence 
of the pair potential as 


A(x, 0) = A(0) tanh (x/^), (A.7) 

with ^ = xf/Aq, the Green functions are represented 

by21 


s(i, e, ^ cosh-“ , (A.8) 

/p(x, 6»,ia;„) =^^tanh , (A.9) 

/i(x, 0,*a;„)=-|^ cosh-2 

where -|- A^(0) and A(0) = Ag cos(0). The 

Green function /p represents the principal pairing corre¬ 
lation in the bulk state, whereas fi represents the pairing 
correlation induced by a surface at x = 0. They are cal- 
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culated from the anomalous Green function as 


fp{x,9,iUJn) 


fi{x,9,iuJn) 


^ (/ + Spf) 

I if- Spf) 


{x,9,iujTi) 

{x,9,iujn) 


(A.ll) 

(A.12) 


At the deep inside of the supercondcutor (i.e., a; :§> ^), 
we obtain fp = f and // = 0. 


1. Current density 

/^From an expression of electric current in Eq. (II3 , we 
define a linear response function TZp^i, by 

jpM = - (A.13) 

= 47rT I ^kpk,yd^^g{r,9,iuJn). (A.14) 

with k = (cos 9, sin 9) and Ue = VpNom being an electron 
density in two-dimension. The diagonal elements of the 
response function correspond to so called pair den¬ 
sity. In the present situation, by substituting Eq. (IA.8I) 
into Eq. (IA.14I) . we obtain 

= 1 — Ki— cosh“^ (, (A.15) 

ne Wo V ? / 

/ rift 9 

-sin2(0)|cos(0)| = -, (A.16) 

where ojq = nT is a low energy cut-off in the Matubara 
summation. Using the normalization condition, the inte¬ 
grand in Eq. (IA.14|) can be represented in an alternative 
way. 


response to an external magnetic field. Eq. (IA.15I) sug¬ 
gests that the paramagnetic response is stronger in lower 
temperature. 

2. Free-energy density 


Substituting Eqs. (jA.8|) - (IA.10l) into Eqs. (I^ - (I2T)) . we 
find that a free-energy density at A = 0 


Tf =NoK2Allog tanh^ , 


Tg = - Non^Al log tanh^ 


K 2 


+ A’oAgK 2 




As a result, we obtain 


Ta 


NoAqK2 



(A.18) 


(A.19) 

(A.20) 


(A.21) 


The free-energy density becomes positive at a; = 0 due 
to appearance of odd-frequency pairs. 

Contribution of a magnetic field to the free-energy can 
be evaluated by applying the expansion in Eqs. (IA.5I1 - 
(IA.6I) onto Eqs. (IA.8I) - (IA.10I) . Within the second order 
expansion, we find that both the pair potential obtained 
by Eq. (jH) and Tf remain unchanged. The second order 
correction to Jy, is given by 


= 2^0 [evFA{x)Y 


1 ^0 1-2 
1 — KI — cosh 
Wo 


(A.22) 


= [-fpdu;„fp + /g. (A.17) 

It is possible to confirm that —fpdai^fpfg corresponds to 
the first term in Eq. (IA.15I) . whereas f\du]^f\/ g contribute 
to the second term. In this way, we can confirm that in¬ 
duced odd-frequency Cooper pairs indicate paramagnetic 


f 2) 

The results indicate that Tg can be locally negative 
(paramagnetic) at low enough temperature. This ex¬ 
plains the decrease of the free-energy in a magnetic field 
shown in Fig. [TOl We also obtain the electric current in 
Eqs. (IA.13I) and (IA.15I) by 2 y{x) = —dTjdAy. 
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